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PROBLEM 46 


Change of Base 


In previous "trip" Problems (The Pi Dragon, The 
Road to e, The Web of Fibonacci), the computing part of @ 
the problem consisted of the moves to be made on the 
coordinate grid. 
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In this problem, the moves themselves are simple, 
and are completely defined by Table A. To insure clarity, 
Table B shows the result of the transformations of Table A 
all applied to the point (-9, -9). 
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Table C (page 4) is a table of factorials, in base 
9 notation. The circled numbers are the low order, 
non-zero digits. The sequence of those digits dictates 
the moves for this trip. That is, we are to apply 
transformations 1, 2, 6, 6, 3, and so on, to X and Y 
coordinates, starting at the origin. Thus we have: 


Transformation 
number 


START 


“N oN ne) Cos (om toy, he) ae 


The first ten moves of the trip are shown graphically 
on the cover. The Problem is this: Where is the 200th 
point? 


The real problem, then, is the calculation of the low 
order non-zero digits of the factorials in base 9. Table C 
shows the first 20 of these. Table D shows the digits from 
101 through 200. The digits from 21 through 100 are to be 
calculated, and all 200 digits used to dictate the transform= 
ations for the trip. 


D Low order non-zero digits of the factorials, base 9; digits 
for positions for 101! through 200! 


43367 45663 47833 66637 15336 33685 76632 24336 74588 73361 
72775 66325 48873 36172 22433 67451 12663 82733 68576 63336 
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fig. C 


Back issues are still available. 
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Given four blocks of storage 


Ten words addressed at 
Ten words addressed at 
Ten words addressed at 
Thirty words addressed 


A Merging 


Problem - 


PROBLEM 47 


as follows: 


A through A+9. 
B through Bt9. 
C through C+9. 
at D through Deg. 


G-t10d 


Blocks A, B, and C contain numbers which are in ascending 
order within each block; there are no duplicates among 
these 30 numbers. We want to merge the 30 numbers into 
block D. (It would be feasible to simply move all 30 
numbers into block D and then sort block D, but this would 
be inefficient. ) 


This is to be a subroutine. The main routine has already 
verified that blocks A, B, and C are as stated, so the 
subroutine need not edit the data. 


A) Draw a flowchart of the logic involved. 


B) Outline a procedure to test a debugged program that 
follows the logic of that flowchart. 


There is shown below a set of sample data, to insure that 
the situation is clear, but of course the logic must apply 
to any data that fits the given conditions. 
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Roots to Order 


In the study of numerical methods, it is expedient 
to have equations at hand for which various algorithms 
for finding roots may be applied. What is wanted are 
polynomials of not too high degree, with integral 
coefficients that are fairly small, and having irrational 
roots. For pedagogical reasons, the roots should be 
easy to predict by the instructor. 


Consider the possibilities: 


1. Use quadratics; the roots can be checked by 
formula. The roots can also be found by formula, and 
a student may properly wonder what we're doing. The 
degree is too small. 


2. Fabricate an equation by building up linear 
factors, such as: 


(oe = (a Fal h/2 ie = Se Os 


Such an equation will either have rational roots, or, 
if it has irrational roots, then its coefficients will 
also be irrational. 
3. Combine the first two methods, as in: 
(x2 4 4x - 1) tx 23) e=E0 


where the left factor has zeros at (-2 +V11). But the 
third root, again, is rational, which spoils the problem. 


4, Use stock equations for which the roots have 
been calculated, such as Wallis! equation: 


x3 o 2s o 5S ©) 
for which the real root is known to some 2000 digits. 

5. Make up an equation with known rational roots 
and then translate it vertically so that the roots become 
irrational. For example, the equation: 

x3 - 5x° = 29x + 105 = 0 
has roots at 3, -5, and 7. Thus, the equation 
x3 - 5x2 - 29x + 104 = 0 


——_———_-- 


has roots that are near 3, -5, and 7 and are irrational. 
(The smallest root is 2.9689.) 


6. Work the problem backwards. Cardan's formulas 
solve the cubic analytically, so we can work from the 
inside out. The critical part of the formulas calls for 


the value of: 
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so we can pick values of p and q to make that term H 

rational. For example, if q is 6 and p is 9, the radical ! 

has the value 6. Then, for | 
As (-q/2) + R 


(-q/2) 


' 
(3.5) 


B 


i 


tt 
‘S) 
= 
m 


a root of the cubic x3 + px + q 
WA + WB 


For the example given, we have: 


W3 + W-9 


which can be readily calculated from a table of cube 
roots: 


- .6378342527444957 322 
for the equation x3 + 9x +620, 


The equation can be translated by replacing x by 
x+k, For example, if x-3 replaces x, we have: 


Xone Qx* + 36x ~ 48 = 0 


for which the roots are 3 greater than the original, or 


X = 2.362165747255504. 
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Gauss's Lattice Problem 


PROBLEM 48 


Tis problem is expounded in a booklet, "Lattice 
Points in a Circle; Experiments and Conjecture," by 
M. E. Rose, Computing and Mathematics Curriculum Project, 
University of Denver, Department of Mathematics, Denver, 
Colorado 80210. 


The problem is this: how many points of a lattice are 
in or on a circle of radius R centered at the origin? 
The Figures show the cases for R= 1, 2, 3, and 4, for 
which the count of points in or on each circle is 5, 13, 
29, and 49 respectively. It was a conjecture of Gauss! 
that it is not possible to write a formula for the number 
of points, Q, as a function of R. For a given R, the 
number Q can be counted by finding all values of X and Y 
that satisfy 


2 + ¥° < Re (1) 


xX 
Clearly, a direct evaluation of (1) would be 
inefficient, since it does not capitalize on the symmetry 
of the problem. Nevertheless, the flowchart (2) shows 
a straightforward approach to the problem. 


Much greater computational efficiencies can be 
obtained by observing (as Rose's paper does) that the 
decision involved in (1) is trivial for most of the points 
within the circle. The points for which computation is 
needed are those lying close to the circle. By using 
this idea and other shortcuts, Richard Sandin calculated 
(1/14/72) the results shown in the following table: 


900 2544569 
1000 3141549 
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radius of circle 
Y: coordinates of 
a point 
intermediate point 
couriter 
number of points in 
or on circle R 


Miriam A. Dumer 
14 January 1972 


(2R+1)°-4N —>Q 


N+1 —>N 
X-1 —>x 


A possible solution to the lattice problem of Gauss. 
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Formulas 


The sum of the consecutive integers from 1 to K 
is given by: 


K(K + iy 


2 


and for the consecutive integers from L to K by: 
(ipo ek > Ee 


The sum of the squares of the consecutive integers from 
1 to K is given by: 


(1/6)(K)(K + 1)(2K + 1). 


The sum of the cubes of the consecutive integers from 
1 to k is given by: 


(2/4)(6")(K +1). 


The sum of the 4th powers of the consecutive integers 
1 to K is given by: 


(1/30) (6K? 4 sk" + 10K° 


- K). 


The sum of the 5th powers of the consecutive integers 
1 to K is given by: 


(1/12)(2K° + 6x? + 5K? - K? 


+ 5K -K) 


The sum of the 6th powers of the consecutive integers 
1 to K is given by: 


6 5 3 


(1/42) (6K! + 21K” + 21K” - 7K~ + kK). 


For the series 


1-3 + 2-4 + 3-5 + 4-6 + 5+7 +...4K(K + 2) 


the sum is given by: 
(1/6) (2x3 + gK° + 7K). 


A Way to Sort 


PROBLEM 43 


The 50 words on the left are 


in alphabetic order. 


is not the 
major sort 


The ordering 
customary one: the 
is on the last letter, 


in descending order; the inter- 
mediate sort is on the second last 


letter, in 
minor sort 
letter, in 


(A) 
accept any 
words and output them with the 
same sorting scheme. 


(B) 


ascending order; the 
is on the third last 
descending order. 


Write a program to 
number of 5-letter 


SATISFY 


The same scheme has 


been applied to the 25 words on 


the right. 


Modify your program 


from (A) to accept any number of 
words of any length greater than 
two letters and output them 
resorted the same way. 


DECISION 
NATION 
ALGORITHM 
MUSEUM 
PRINCIPAL 
ASTONISH 
STRING 
SHERIFF 
TIDE 
STORAGE 
STRIVE 
TOLD 
BASIC 
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Sequences of Triangles 


PROBLEM 38 SOLUTION 

Starting with an equilateral triangle with unit 
area, a sequence of equilateral triangles is formed in 
which the area of one is the altitude of the previous 
one. The problem was to find the altitude of the 100th 
such triangle. 


Several readers pointed out that the sequence 
converges, so that 


=m «4 / 
en =V3 Ay 
from which it can be deduced that 
i. SW s3ho 


ioe) 


But that wasn't the problem; the problem was to find Aigo: 
Associate Editor David Babcock calculated: 


A109 = 1+73205080756887729352744634 15051218 


which agrees with the square root of 3 (see PC3-6) to 
31 significant digits. 
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Solution 


Problem 11 (PC5-13), the 
100-Square Trip Problem, 
called for the logic of 
traversing a 10 x 10 
array up and down its 
columns. The subroutine 
logic on this flowchart 
was devised by Tom Cundey. 
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THE | 


|Maze Game 


PROBLEM 50 


Given the numbers from 1 to 100. Starting at l, 
a move is made to another number according to these rules: 


From number X (1) 
proceed to (2) 
the number (3) 
given by: (4) 
(5) |x? - x3 


[all arithmetic is modulo 100, and rule (5) operates on the 
results of modulo 100 arithmetic for the squares and cubes] 


One of the five rules will usually select a new number; that 
is, one not previously chosen. If not, then the lowest 
number still available is selected. When the number 100 

is reached, the game is over. Normally, a11 100 numbers 
will be selected in each game. 


Consider the selection of the numbers as a journey 
on the pattern shown here: 


73| 74 | 75|76 [77 | 78 | 79 | 80) ¢ 

ra | | | |an 
42 [2 | 22 | 25 | 24 | 25 |26 [81 [84 
j4i}20] 7 | a | 9 [10 |27 [52 | 25) 

[40 fag] 6 | 1 | 2 | 11] 28 §3 | a6 


6a {sa]te] 5 [4 | 3 [12 |29 [54] 7 | 
67 | 30 [17 [16 [15 [14 [13 | 30] 55 | 
66 | 37 | 36 | 36 | 34 | 33 | 32 | 31 | 56/69 
65 64 | 63 | 62] 61 | 60 | 5a | 88 5790 | 
t00| 99 88 | 97] 96 | 95 | 94 93) 92) 91 | 


the distance traveled can be calculated. For convenience, 
the square of the distance is used. Thus, the D“ distance 
from cell 1 to cell 84 is 5° + 22 = 29, ‘The journey 


begins at 00, so the first leg, to cell 1, has a D¢ 
distance of 50 to start. 


The Problem is, what ordering (following the move 
rules) will produce (A) the longest journey, or (B) the 
shortest journey? Present records are 2711 for A and 
1803 for B. The shortest known journey begins as 
follows: 


arbitrary 
only possible move 
x3 


Gaifference between powers 
by rule (6)--by default 


Psa 
difference 


Solution 


Problem 37 (PC12-1), the Sine Excursion trip, 
called for a 600-leg journey in which the lengths of each 
leg were given by the decimal expansion of sine 1, and 
the turns were uniformly one radian clockwise. 


Thomas R. Parkin, Control Data Corporation, furnishes 
these results: 


x 32 .9678624079 
Y = -70.4643240552 
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N-Series 


1. 1461280356782380259259551533171292202517622777860 
739478140624.14845361629 17650367 555303877996567475 


2.639057 329615258614.5225848649013562977125848639421 
164425800701594 3097 348472176398 339 352182558429021 


3.741657 386773941 3855837487 32316549 3017560198077787 
26946 303745467 32003515630693902797 6809895194.37958 


2.410142264175229986128 36966760327 2895 3545812899808 
6765416413971041 329172692259 383382261151622681347 


1.6952182030724 35481549 34 3584607 767115294 3805646840 
91593099616 35805458 32 36090808 1774415890032537 1200 


1.457916249576283530691311271122606934 3069267644713 
5425221119466449 337925197 185565657078460176015252 


1. 302005454317467 704497249 3030774256 3032 30288915111 
935397627 184827 375737757 098509914886787 3589479168 


1 .0267418881337292354684 536 395104159442 321062634164 
5761923285260174114929108109109452 3484414 36523084 


1202604 .2841647767777492 367707 6785944941248654 33761 
0224031329063319746294708334267090364192964 


9122171.1817543531702043751107628162745027008832977 
6225299 3768387 3097427 636237 77951986 3008 3460 


1.4994888620096062927989507017866583810752847684575 
108316742798 32024 36565297817 683027845 30268807 1088 


41001860888499 3288052964. 1652467 0972545801067 5237920 
27 322197 126356748926146602648 3061479032219018658198 
1413953765376 


